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Abstract 

We propose a nonlinear integral equation (NLIE) with only one 
unknown function, which gives the free energy of the integrable one 
dimensional Heisenberg model with arbitrary spin. In deriving the 
NLIE, the quantum Jacobi-Trudi and Giambelli formula (Bazhanov- 
Reshetikhin formula), which gives the solution of the T-system, plays 
an important role. In addition, we also calculate the high temperature 
expansion of the specific heat and the magnetic susceptibility. 
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1 Introduction 

Thermodynamic Bethe ansatz (TBA) equations have been used to analyze 
thermodynamics of various kind of solvable lattice models p. However it 
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is not always easy to treat TBA equations since in general they contain an 
infinite number of unknown functions. So there are several attempts to sim- 
plify the TBA equations. In particular for the ID spin | Heisenberg model, 
Kliimper proposed nonlinear integral equations (NLIE) [2j which contain a 
finite number of unknown functions from the point of view of the quantum 
transfer matrix (QTM) approach U 13 121 • There is also a similar 
NLIE by Destri and de Vega |H]. 

Another type of NLIE for the spin | XX Z model, which contains only 
one unknown function, was discovered by Takahashi ^ to simplify the TBA 
equation. Later, this NLIE was rederived ^2] from the T-system [TT| IT^ of 
the QTM. Moreover this NLIE for the spin ^ XXX model was also rederived 
from a fugacity expansion formula JHl- We derived this type of NLIE for 
the osp{l\2s) model p-4] and the sl{r + 1) Uimin-Sutherland model J3] for 
arbitrary rank. The number of the unknown functions of these NLIE coin- 
cides with the rank of the underlying algebras. All of them are NLIE for 
fundamental representations of underling algebras. So it is desirable to de- 
rive NLIE for tensor representations, i.e. NLIE for higher spin models. The 
purpose of this paper is to derive NLIE for the Heisenberg spin chain with 
arbitrary spin | [16 . 

Thermodynamics of the higher spin Heisenberg model was firstly inves- 
tigated [16j by the TBA equations, which contain an infinite number of un- 
known functions. Later a set of NLIE with s -|- 1 unknown functions for this 
model was derived ^Tj by the QTM approach. This NLIE is an extension of 
the Kliimper's type of NLIE [2^. On the other hand, our new NLIE contains 
only one unknown function for arbitrary s. Thus, as far as the number of 
the unknown functions, our new NLIE is a further simplification of the TBA 
equations. 

In section 2, we introduce the higher spin Heisenberg model and define 
the T-system of the QTM. As a solution of the T-system, we introduce the 
quantum Jacobi-Trudi and Giambelli formula (Bazhanov-Reshetikhin for- 
mula 18 J ()2.18|1 . which plays an essential role in the derivation of the NLIE. 
This formula expresses an eigenvalue formula of the transfer matrix for the 
tensor representation in the auxiliary space in a determinant form. In the 
representation theoretical context, it may be viewed as a Yangian analogue of 
classical Jacobi-Trudi and Giambelli formula on Schur functions fHll^- In 
section 3, we derive the NLIE ()3.9|1 . which is our main result. The T-system 
which we have to use here is not the standard one ()2.15|) but an old one 
[TlJ (|2.20|) . ()2.21|) . Due to the quantum Jacobi-Trudi and Giambelli formula, 
determinants appear in our new NLIE. These novel situations contrast with 
the fundamental representation cases [HI EH El E|- Using our new NLIE 
()3.9j) . we also calculate the high temperature expansion of the specific heat 
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and the magnetic susceptibility in section 4. It will be difficult to obtain 
the same result by the traditional TBA equations. Section 5 is devoted to 
concluding remarks. 



2 QTM method, T-system, quantum Jacobi- 
Trudi and Giambelli formula 

We introduce the higher spin Heisenberg model, and define the QTM, the T- 
system and the quantum Jacobi-Trudi and Giambelli formula for this model. 
A more detailed explanation of the QTM analysis for this model can be found 
in Hlj. 

The Hamiltonian of the spin | Heisenberg model is given as follows 

L 

^0 = JJ]Q.(S,S^.+i), (2.1) 

i=i 

where Sj = {Sj, Sj, Sj) is the spin | operator acting on the j-th site, and 
SjSj+i = SjSj^i + SjSj^-^ + SjSj_^_i. J is a real coupling constant (J > 
(resp. J < 0) corresponds to the anti-ferromagnetic (resp. ferromagnetic) 
regime) and L is the number of the lattice sites. We assume the periodic 
boundary condition Sl+i = Si. Qs{x) is defined as 

s.w = -2EEi n (2-2) 

j=0 k=j+l p=0{p^j) ^ 

where Xp = |(p(p + 1) — s(| + 1)). For s = 1 and s = 2 case, ()2.1|1 becomes 

L 1 J ^ 

Ho = 2J ^{SjSj+i - -}, ^0 = ^{SjSj+i - (SjSj+i)^}, (2.3) 

j=i i=i 

respectively. The i?-matrix [201 1^ of the classical counter part of the spin 
I Heisenberg model is defined as 

^i^") = E( n "i^)p'' (2-'') 

j=0 k=j+l 

where is the projector onto j + 1 dimensional irreducible s/(2) module, 
can be expressed as 

P== ff ^^^-'^v ^ ^2.5) 



where S = {S^, S^, S^) is the spin | operator, and S ® S = ® + ^ 
S'^ + ^ . In our normahzation of the i?-matrix, R[oo) is an identity 
operator. The row to row transfer matrix is defined as 

t{v)=troRoLiv)---Ro2iv)Roiiv), (2.6) 

where Roj{v) is the i?-matrix ()2.4j) acting on the auxihary space and the 
j-th site of the quantum space. This transfer matrix is connected to the 
hamiltonian ()2.1|) as 

no = J^\ogt{v)U=o. (2.7) 

One can add an external field term Tiex = —'^h Sj in the hamiltonian 

1)2.11) without breaking the integrability as = TYq + 'Hex- The QTM is 
defined as 

2hSg ~ ~ 

tQTM{v)=troe T Ri^q{u - iv)RN-i o{u + iv) ■ ■ ■ R2o{u - iv)Rio{u + iv){2.8) 

where Rjkiy) is defined by '90° rotation' of Rkj{v), i.e., Rjkiy) = ^'^Rkjiv) 
{tk- the transposition of Rkjiy) in the k-th. space); is the Trotter number 
and assumed to be even; u = (T is a temperature); the Boltzmann 

constant is set to 1. The free energy per site is expressed in terms of the 
largest eigenvalue fs{v) of the QTM at v = 0: 

/ = -T lim logT;(0). (2.9) 

We can embed Ts{v) into the eigenvalue formulae {T!m(f )}mez>o for the fusion 
hierarchy of the QTM, i.e. a Yang-Baxterization of the character of the m-th 
symmetric tensor representation of s/(2): 

f^{v) =T^{v)/Arm{v), (2.10) 

where 

m 

Tm{v)= Yl l[z{df,v+'-{m-2j + l)) 

l<di<d2<---<dm<2 3=1 
m rri-k „ . o ■ i o 

E (m-fc)Mi+fcM2 -TT , , m — s — 2] . ^ , m ~ s ~ 2] + 2 . 
e T y[ + -t)M^ + 77-^ ^) 



k=0 j=l 

h 
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n, , m + s - 2j ., , , m + s - 2j + 2 , 
0-(^ + ^ — -^)M^ + ^-^ — ^) 

j=m—k+l 
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. r / N -TT / m — 2A; — s,,, m — 2k + s + 2 ^ 

M^{v) = \[ct>^{v + + ^ 0; (2.12) 

fc=l 

0±('i^) = {v^iu)^] Q{v) = nj=i('^~'^i) ^ ^>o); A^i and are chemical 
potentials (in our case, /ii = h, fi2 = —h); 



z{\\ v) = e T — — i 



2 ^'g(t;-f)' 

z{2,v) = e^0_(- + + • (2.13) 

Here we adopt a convention nj'=i(' ' ' ) = njlm+i(' ' ' ) = L "^^j is a solution 
of the Bethe ansatz equation 

.#._(t.t + f)M"t + f + ¥) + i) 

where g = 2: the dual Coxeter number. Note that T^iv) ()2.1ip is free of 
poles under the Bethe ansatz equation ()2.14|) . One can show that the function 
Tm{v) ()2.1()|1 satisfies the following T-system [T2j. 

Tm{v - ^)Tm{v + ^) = Tm-i{v)f^+i{'") + 9m{v) , m G Z>i, (2.15) 

where 

To{v) = 1, (2.16) 

min(m,s) . m+s+l-2fc - / _ m+s+l-2fc ■^ 

The solution of the T-system ()2.15p is given by the following quantum Jacobi- 
Trudi and Giambelli formula (Bazhanov-Reshetikhin formula |18j). 

where the matrix elements are given as follows 

1 a = 

f'^^'^ = < It) T=l (2,19) 

a > 2 or a < 0. 



The following functional relations ^T] are equivalent to the T-system ()2.15|) . 

~ ~ m ~ m — 1 _ i ~ m + 1 
Ti[v)Tm-i[v - —%) = Tm[v ^ — i) + gi{v - -)Tm~2[v ^ — (2.20) 

T^{v)T^^x{v + -%) = Tm{v + ^) + 9i{v + -)T™_2(t^ + 0> (2-21) 

m G Z>i. 

()2.20|) follows from an expansion of the determinant ()2.18|) down the first 
row. p.21|l follows from an expansion of the determinant ()2.18|1 down the 
last column. 



3 A new nonlinear integral equation 

Suzuki had J7j the following observation from numerics. 

For the largest eigenvalue sector ofTs{v), the corresponding root of the 
BAE \2.14^ forms s-strings. In this case, imaginary parts of the zeros {5^} 
ofTm{v) (m = 1, 2, . . . , are located near the lines Imv = ±|(s + m — 2j), 
j = 0, 1, . . . ,m - 1. 

Following Suzuki's calculations for the spin | = 1 case, we have plotted 
roots of the BAE ()2.14j) in the sector M = N and zeros of Tm{v) (m = 1, 2) 
(FigureC]and FigureE))- Admitting Suzuki's observation, we shall derive the 
NLIE for the largest eigenvalue sector of Ts{v) from now on. Tm{v) has poles 
1 at = ±/3^ : = + u)t, + u)t, ( "^+^+^-2^^^^^) ^ ^^^^ 

whose order is at most ^. Moreover 



m m 

{m— fc)/j.]^+fc/i2 — ^ h{m — 2k) 



:= hm TUv) = V e = V (3.1) 

' ' k=0 k=0 



is a finite number. This is a solution of the Q-system 

(Q™)' = Qm~iQm+i + Q'^\ m e z>i, (3.2) 

where Qm'* = e ^''t = 1 and Qo = 1. So we must put 

— _ 

Ti{v) = Qi + J2 {-At- + -At- \ ^ (3-3) 

7 = 1 



^Note that these poles are known ones, which we need not solve the BAE H2.14|l to get. 
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Figure 1: Location of the roots of the BAE ()2.14|) for the spin f = 1 case 
(A^ = 12, u = —0.05, J = 1, /i = 0, 1, 2), which will give the largest eigenvalue 
for T2{v) at f = 0. The root forms 2-strings at least for h = 0. 
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(m,h) = (1,0) o 
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Figure 2: Location of the zeros of Tm{v) (m = 1,2) for the roots in Figure 
HJ The zeros are located outside of the physical strip Imt> G [— |, |]- 
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where jSi = {^-^ + u)i. Here the coefficients are given as follows. 

h, = j^^Uv){v - h = £^T,{v){v + P,y-\ (3.4) 

where the contour C (resp. C) is a counterclockwise closed loop around (3i 
(resp. —Pi) that does not surround — /3i (resp. Pi). Substituting ()3.4|) into 
(nO|l . and using (jT^ and ^TH^ . we obtain 



Ti(t;)=gi 



^ 9i{y+ ~Pij i)Tm-2{y + A - 
fm-i{y + Pi-ft) 

JV 



c27r^ v-y + Pi tT„_i(y-/?i + fz) 

^ 5^i(y-/3i^+f)T^-2(y-/3i + ^^) l 
T^_i(y-A + fO i' 

where the contour C (resp. C) is a counterclockwise closed loop around 
that does not surround —2Pi (resp. 2Pi). The ffist term and the second 
term in the first bracket {■ ■ ■ } in ()3.5p have a common pole at 0. However 
this common pole of the ffist term disappears if m > s + 1. Therefore, for 
m > s + 1 the ffist term vanishes after the integration as long as the contour 
C does not surround the pole at Zm-i — Pi + ^i- Similarly, for m > s + 1, 
the first term in the second bracket { ■ ■ • } vanishes after the integration as 
long as the contour C does not surround the pole at z^-i + Pi — ^i- Thus, 
for m > s + 1, we obtain 



Ti{v) = Qi + 



dy gi{y + - t)T^-2(y + A - ^i) 
elm {v-y- ~Pi)Tm-i{y + A - f 



2 

^ . dy ^i(y-/3i + f)T^-2(y-/3i + ^z) ^gg^ 



c 



2m ^y-y + p^)T^_^^y-Pi + ft 

JV 

Here we omit the terms which contain y~ since these terms cancel the poles 
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of Qi. We can take the Trotter limit N oo. 

dy gi{y + A - f )T^^2(y + A - ^i) 



^2m {v-y- Pi)r^-i{y + /3i - f ^) 
dy gi{y - A + f )T^-2(i/ - A + 



(3.7) 



% 2712 iv-y + (3i)%n.iiy - A + f z) 
where (3i = hniAr^ooA = %n{v) = limN^ooTm{v) and 

^iM=exp + =Qf^exp l^—^^l. (3.8) 

In particular, ()3.7|1 for m = s + 1 is the simplest. 

d?/ gi{y + ^i)rs^i{y - ^) 



c27rz (v - y - ^i)%iy) 



dy 9i{y - ^i)Z-i{y + ^) 
%27n {v-y + ^t)Z{y) 



where the contour C (resp. C) is a counterclockwise closed loop around 
that does not surround —(1 + s)i (resp. (1 + s)i) and Zg. ()3.9j) contains 
only one unknown function Ti{y) since Tsiy) and 7^_i(f) can be expressed 
by Ti{v) through ()2.18j) in the Trotter limit. The free energy per site is given 
by / = -riogT,(0). For s = 1, dSH) reduces to the Takahashi's NLIE for 
the XXX spin chain Although we only consider the largest eigenvalue 
of Tsiv) in the limit — oo , other eigenvalues also satisfy the NLIE ()3.9|) 
if above conditions for the integral contours are satisfied. In usual, such 
eigenvalues have zeros in the physical strip \mv G [— |, |]- Thus to exclude 
the eigenvalues other than the largest one, one may take the integral contours 
on the line Imv = 



4 High temperature expansion 

In this section, we shall calculate the high temperature expansion of the free 
energy from our new NLIE ()3.9|1 . The calculation is not easier than s = 1 
case j21j due to the determinants in ()3.9|1 . However it is easier than to use 
the traditional TEA equation. For small J/T, we shall put 

r^(t;)=exp(f;6„,„(t;)(^)"), (4.1) 

\n=0 / 
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where 6m,o(^) = logQm- Due to ()2.18p in the hmit ^ oo, one can express 
bm,n{v) in terms of fundamental ones {bi^k{v)}o<k<n, Qf^ and h{v) = s/(t>^ + 
For example, we have 



b2Av) = [ iQi)X,iiv -'-) + iQi)%,iiv + J) - Q?biv) ] /Q2, (4.2) 



where Q2 = [QiY — Qf^- Taking note on the relations like ()4.2|1 . substitute 
()4.H) into we obtain the coefficients {bm,n{v)}. For example, {bi^n{v)} 

for s = 2, n = 1, 2, 3 are as follows. 



bi,i{v) - 



+ ((gi)2-gf)^ 

b,,2iv) = 2gS')| (45 + Av^) (Qi)^ + (-99 + Av^) (QiYQ? 

-4 (27 + 20.^) iQf^r}/ I (9 + Av')\iQ,r - Q?)'^ , 

biM = {Qf^ ((477 + 120^;' + 16v^) (Qi)« + 9 (-179 + 5Qv^ + 16v^) {QifQ^^^ 
-6 (1101 + 1536v^ + 30Av^) {Qi)\Q'-^^y 

+4 (5553 + 5028v^ + 89Qv^) {Qi)\Q?f + 48 (63 + 84i;^ + 32i;^) (Qf ^)^) } 
/{{9 + ^v^)\{Qir-Q?y]. (4.3) 



We can calculate the specific heat C = —T-^^ and the magnetic suscepti- 

bility X = —-Qf^\h=o- In this case, we only use {bm,k{v)} for t> = due to 
the definition of the free energy ()2.9j) . Note that the /i-dependence of the 
free energy enters only through Qi = +e~^ since Q^^ = 1. Let us put 
t = J/T. 
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s = 2 case: 

^ Ut^ 580^5 27629 165529 40875277 1« 
C = — + — — -— — + + 



9 27 54 243 11664 116640 8398080 
1064887^^^ 2176205977 94355582827 100181936647 



4408992 470292480 12697896960 304749527040 

23030455724107 4273471680238097 1^'^ 919571868309546869 

^ 2394460569600 ^ 482723250831360 188262067824230400 
30226224111769481353^16 1741566568691332074437 



1916850145119436800 237210205458530304000 
636000568233936625686389 1^^ 5325661218179817974957 1^^ 

+ 45544359448037818368000 ^ 248158368787385548800 
143978522302549633307591^20 3868537325098539118572347677^21 



^ 227535901895503223193600 144433584998668295995392000 
+0(t''), (4.4) 

_ 8 8t lAt^ 205 97 32627 1"^ 290839 2993083 

yT = \ \ \ \ 

3 3 27 324 810 58320 489888 26127360 

51476713^9 392473169^1° 6147244063 28552400626009 

^ 70543872 ^ 846526464 14549673600 33522447974400 

41 136528139217^13 46001003925515^1^ 5761065533176715581 f^^ 



186767924428800 59765735817216 6589172373848064000 
393055654682062161341 t^^ 1 142327068285920573167 



2530242191557656576000 1024145648963813376000 
+0(ti«). (4.5) 

s = 3 case: 

^ 15^2 145^3 1385 13445 1203755 7458199 
C = —— + —— + 



16 96 2304 6912 331776 5971968 

66731599 29819652545 548487273433^1° 62192420939825 

^ 15925248 ^ 4514807808 ^ 433421549568 7801587892224 
29359344590299711 t^^ 1535815314115199ti3 



2808571641200640 4413469721886720 
12355531497035499295^1^ 5343429626131816107323 nr lew. 
^ 889755495932362752 ^ 347004643413621473280 ^ ^'^ ' 

35^3 2725 3383 173425 7164095 
xT = 5 - 5t +—— + —— + 



54 2592 5184 373248 5878656 

9795049 i« 634647695 2120381481515 13785996387509*11 
+ ^.^^.^^^^ + , ^ ^ + 



13934592 952342272 1462797729792 20113468784640 

33615348883175267*12 37654994203903963 *i3 ^,,4, 

+ 0(^1^). (4.7) 



34756074059857920 21515664894197760 
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s = A case: 



24^2 619 49159 900239 257362861 

~ ^5~ ^ 375 ^ 45000 675000 64800000 

190318307851 1061704692647 33531890711924393 

+ + 



+ 



58320000000 518400000000 4408992000000000 
13698973673330960069 1^^ 4091902458911705383 



2116316160000000000 1360488960000000000 
1767044044952349905869283 1^^ 



137137287168000000000000 



+ 0(t"), {U 



101 9337 209797 7464779 



xT = 8 - 8t +—— + —— + 



135 6480 162000 116640000 
234988285877 3029100708947 393143129330767^^ 



146966400000 1741824000000 3809369088000000 

4303789688178371831 278589255303992797247tii i2x.. ^ 

+ 0(t^2).(4.9) 



2285621452800000000 125709179904000000000 

We have plotted ^ the high temperature expansion of the specific heat and 
the magnetic susceptibihty in Figure 01- Figure El According to the Figure 
El the position of the peak of the specific heat seems not change drastically 
when s changes. In particular, s = 2, 3 cases agree with the result from 
another NLIE for large T (see Figure 6 in PZ]). This indicates the validity 
of our new NLIE. 



5 Concluding remarks 

In this paper, we have derived a NLIE with only one unknown function. 
This type of the NLIE for higher spin Heisenberg model with arbitrary spin 
is derived for the first time. In particular, a remarkable connection between 
the NLIE and the quantum Jacobi-Trudi and Giambelli formula is firstly 
found. 

Now that the NLIE is given, the next important task is to search the 
solutions of it. The zero-th order of the high temperature expansion 1)4.111 

^Here we have used the Pade approximation. The [L, M] Pade approximation of a func- 
tion g{t) of t is the ratio of a polynomial of degree L and M: g{t) = + 
Q^^L+M+i-j^ It provides approximately an analytically continued function of g{t) for out- 
side of the radius of convergence of g{t) . Thus the Pade approximation is expected to 
provide better results for small T than original plain series. For more detail, see for 
example, |2S1- 
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Figure 3: Temperature dependence of the high temperature expansion of 
the specific heat C for J = 1 at /i = 0: broken hues denote plain series 
and smooth hues denote their Fade approximation for s = 1, 2, 3, 4 from the 
bottom to the top. Their orders are (n: plain series 0(1/T"'), Fade) = (52, 
[26,26]), (21, [10,10]), (15, [7,7]), (12, [6,6]) respectively. Their peak positions 
and peaks are (peak position, peak) = (0.961, 0.350), (0.963, 0.589), (0.956, 
0.780), (0.949, 0.942) respectively. The case for s = 1 was calculated in ref. 

m 
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T 



Figure 4: Temperature dependence of the high temperature expansion of the 
specific heat C for J = 1, s = 2: broken hnes denote plain series and smooth 
hues denote their Pade approximation for /i = 0, 2, 4 from the bottom to the 
top on the right side. Their order is (n: plain scries 0(1/T"), Pade) = (17, 
[8,8]). Their peak positions and peaks arc (peak position, peak)= (0.964, 
0.589), (1.365, 0.520), (3.452, 0.612) respectively. 
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2 . 5 



1.5 



. 5 




Figure 5: Temperature dependence of the high temperature expansion of the 
magnetic susceptibihty x fo^' </ = 1 at /i = 0: broken hues denote plain series 
and smooth hues denote their Fade approximation for s = 1, 2, 3, 4 from the 
bottom to the top on the right side. Their orders are (n: plain series 0(1/T"), 
Fade) = (28, [13,13]), (18, [8,8]), (14, [7,7]), (12, [6,6]) respectively. Their first 
peak positions from the right and the peaks are (peak position, peak) = (1.282, 
0.294), (1.454, 0.748), (1.555, 1.367), (1.624, 2.153) respectively. The case 
for s = 1 was calculated in ref. 
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is a solution of the Q-system. Thus this task is to incorporate the spectral 
parameter into a solution of the Q-system, namely to find a solution of the 
T-system. The solution of the Q-system is a kind of a generalization of the 
hypergeometric function (cf. |^). Thus we expect that the final answer will 
be a further generalization of the hypergeometric function. If a hypergeo- 
metric series solution is found, one should consider an integral representation 
of it; then one may be able to treat the low temperature region where the 
plain series does not converge by an analytic continuation. 

Finally, we note that we can also derive a NLIE similar to ()3.9|) for the 
row to row transfer matrix. 
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